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Abstract 

The inverse problem of the calculus of variations asks whether a given 
system of partial differential equations (PDEs) admits a variational for- 
mulation. We show that the existence of a presymplectic form in the 
variational bicomplex, when horizontally closed on solutions, allows us 
to construct a variational formulation for a subsystem of the given PDE. 
No constraints on the differential order or number of dependent or inde- 
pendent variables are assumed. The proof follows a recent observation of 
Bridges, Hydon and Lawson and generalizes an older result of Henneaux 
from ordinary differential equations (ODEs) to PDEs. Uniqueness of the 
variational formulation is also discussed. 

1 Introduction 

Many systems of partial diflterential equations (PDEs for short) that appear in 
physics are variational. That is, they are equivalent to Euler-Lagrange (EL) 
equations of some local Lagrangian density. The question of whether some 
arbitrarily given PDE is variational is known as the inverse problem of the 
calculus of variations. It has attracted a significant amount of attention in the 
past [2j[21]. In physics, a variational formulations endows the algebra functions 
on the phase space of a classical theory with Poisson structure, which in turn 
determines the corresponding commutation relations in quantum theory. One 
important consequence of the inverse problem is the ability to reconstruct a 
variational principle and hence the quantum commutation rules (whether unique 
or not) directly from the equations of motion [24[ [15] . 

This problem comes in two versions, the harder multiplier problem and 
the easier non-multiplier problem. The idea behind the solution of the non- 
multiplier problem for ordinary differential equations (ODEs) was known al- 
ready to Helmholtz and the necessary and sufficient conditions for a positive 
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solution are known as the Helmholtz conditions. The generahzation of these con- 
ditions to PDEs is elegantly stated in terms of the variational bicomplex [3l [4] . 
The multiplier problem is much less understood, with significant results ob- 
tained for ODEs by Douglas [12] and for PDEs by Anderson and Duchamp [5]. 
A criterion certifying a positive solution in the case of second order ODEs was 
given later by Henneaux [TS] , which was later extended to higher orders [T71 [T] . 
This criterion is essentially the existence of a (pre)symplectic form on the space 
of dependent variables that is conserved by the flow of the ODE. 

The work of Henneaux proved difficult to generalize to partial differential 
equations (PDEs). For one thing, it is not immediately clear what is the right 
analog of the symplectic form. Also, Henneaux represented a second order ODE 
as a vector field on the tangent bundle of the configuration space. Finally, the 
conservation condition was expressed as annihilation by the Lie derivative of 
the ODE vector field (with suitable generalization in the time dependent case). 
However, we can now say that Henneaux's argument can in fact be (at least 
partially) generalized to the PDE case. The analog of the symplectic form is 
the covariant presymplectic current (Section[3|). And the analog of the geometric 
vector field representation of an ODE is the geometric representation of a PDE 
as a submanifold of the jet bundle of the dependent variables (Section[5]). Within 
this context, the corresponding generalization of Henneaux's criterion certifying 
a positive solution of the multiplier inverse problem for PDEs was given only 
very recently in [SI 116] . This key observation appeared as a side remark in that 
work, with only a sketch of the details and without placing it within the context 
of the inverse problem. Unfortunately, this solution of the inverse problem is 
only partial, as the Lagrangian density produced by the proposed procedure 
(Section [4|) produces a variational system that need not be equivalent to the 
original PDE system, though this variational system will admit all solutions of 
the original PDE system among its own solutions. 

This note aims at highlighting and clarifying the above result by placing it in 
the appropriate geometric context. Section [5] introduces some basic background 
on the geometric formulation of PDEs using jet bundles and the associated 
variational bicomplex. Section [3] shows how a conserved presymplectic current 
density arises in variational problems. Conversely, Section U] shows that the 
presence of a conserved presymplectic current density is equivalent to the fact 
that the solutions of the given PDE system are also solutions (but perhaps 
not the only ones) of a variational PDE system. Section [5] discusses the non- 
uniqueness of the Lagrangian constructed in Section (4] Finally, conclusions are 
in Section [SI 

2 Jet bundles and PDE systems 

This section briefly defines the basic notions needed for the geometric formu- 
lation of PDE systems in terms of jet bundles, including associated geometric 
structures, and fixes some notation. For a detailed discussion of jet bundles 
see [201 HH] ■ The variational bicomplex is discussed in detail in [3l H] . 
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Let F ^ M he a vector bundle over the base space M, an n-dimensional 
smooth manifold. Since we will be only interested in local questions, assume 
that both F and M are topologically trivial (contractible). The arguments 
in this paper can be straightforwardly generalized to smooth bundles with non- 
trivial fiber and base space topologies [T^- Let r(F) denote the space of smooth 
sections of F. Denote by T*N the cotangent bundle for any smooth manifold N 
and by A'^N = /\^j^T*N the bundle of alternating fc-forms over A''. Denote by 
Vi^{N) = T{K^N) the space of differential forms on iV, by d the corresponding 
de Rham differential and by U,*{N) the total complex of differential forms graded 
by degree. 

The fc-jet of a section (j) G r(F) at a; G Af can be defined as the equiv- 
alence class of sections ■0 S r(_F) that have coinciding Taylor polynomials up 
to and including order k at point x in any (bundle adapted) local coordinate 
system. Thus, /c-jets are a coordinate invariant way of capturing the derivatives 
of a section up to order k. Local coordinates (x\ u") on F naturally extend to 
local coordinates (x^wf), \I\ < k on J^F, where / with |/| = Hs a multi-index 
(iiZ2 ■ ■ • ii). The coordinates of the fc-jet of (/) at a; are 

= \I\ = k, (1) 

where we have used the shorthand dj — dt^ ■ ■ ■ di^ . These fc-jets form the space 
J^F, which is a fiber of the vector bundle J'^F M of fc-jets over M . The 
0-jets are identical with the underlying bundle, J'^i^ = F . Jet bundles come 
with natural projections J^F — j'F, for any k > I, which simply discard the 
information about all derivatives higher than I. This projection gives J'^F the 
structure of an afhne bundle over J''F. The projective limit 

J°°F = ^ J'^F > J^F ^ J^F J°F M (2) 

k 

is called the oo-jet bundle. A smooth function on J°°F is the puUback of a 
smooth function on some j'^F for some fc < oo. That is, a smooth function on 
J°° F always depends only on finitely many components of an cxD-jet given to it 
as argument. 

A section (j) G r(F) naturally gives rise to the section ^^(j): M J^F, with 
j''(f>{x) = j^cj), called the fc-jet prolongation of (f>. Similarly, a (not necessarily 
linear) bundle morphism /: J^F E (which determines a differential operator 
f o '■ r(F) — > T{E) of order fc), naturally gives rise to a bundle morphism 
jk+ip _^ jij^^ f{3t^^4>) = jiif ° j''{<l>)), called the Z-prolongation 

of/. 

In the geometric formulation, a PDE system on F of order fc is a submanifold 
l: £ C J'^F that satisfies the regularity conditions of being closed and that 
£ M is a. smooth sub-bundle of J'^F M. To connect with the usual 
notion of a PDE, we note that there always exists (at least up to possible 
global topological obstructions [14, §7]) a vector bundle E M and a bundle 
morphism /: j'^F — > E such that a section (f) G T{F) satisfies the system of 
differential equations foj''{(j)) = iff the image of j''(p G T{J''F) is contained in 
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£, that is j'^cj) is actually a section of £ -t- M. We call the pair (/, E) an equation 
form of £, which is in general not unique. The equation form can always be 
chosen to be regular, which means that any smooth function L on J^F vanishes 
on £ iff it can be written, in local coordinates (x', u") and (x*, v^) on F and E, 
as L{x'-,uf) = LB{x^,uf)f^{x^,u'}), with smooth coefHcients Lb- Conversely, 
any regular equation form (/, E) defines a PDE system £f given by the zero 
set of /. The /-prolongation i/ : C J^^^F is defined as the PDE system £pi j 
corresponding to the equation form (p'/, ,PE). The /-prolongation comes with 
with a natural projection pi: £^ — )• £, which simply restricts the appropriate 
projection of jet bundles. Note that pi need not be surjective if the PDE system 
has non-trivial integrability conditions. Below, we deal with oo-prolongations 
Loo - £°° C J°°F and J°°E)- We assume that both are regular. 

The de Rham differential acting on il*{J°°F) can be naturally written as 
the sum 

d = dh + dv (3) 

of the, respectively, horizontal and vertical differentials. Each is nilpotent and 
they anticommutc: 

d^ = d^ = 0, dhdv + dvdh = 0. (4) 

If d-M is the de Rham differential on M, (j) g r(F) is any section and L G 
^1*{J°°F) any differential form, then the defining property of the horizontal 
differential is that it is intertwined with dM via the puUback of differential 
forms along sections: 

{r^rduL = dM{r^rL- (5) 

Differential forms on J°°F have the following natural subspaces: the purely 
horizontal forms 17'''° (F) that are generated by puUbacks of forms from n^{M) 
along the natural projection J°°F M and the purely vertical forms 17°'" (i^) C 
n'"{J°°F) that are annihilated by the pullback of any section M ^ F- 

Purely horizontal and purely vertical forms generate n*{J°°F) as a graded alge- 
bra. The subspaces of homogeneous horizontal and vertical degrees are denoted 
by n^'^iF) C r2''+^(J°°F). The differentials dh and dv are then, respectively, of 
purely horizontal degree 1 and of purely vertical degree 1. In local coordinates 
(a;*,Uj), horizontal forms are generated by dux^ = da;*, while vertical forms are 
generated by d^Uj- The total bi-differential bi-graded algebra (17*'*(F), dh, dy) 
is called the variational bicomplex of the vector bundle F — M. Within this 
bicomplex, we can define the horizontal, 7?*'*(dh), and vertical, H*'*{dv), coho- 
mology groups in the obvious way. 

The horizontal and vertical degrees, as well as differentials, survive restric- 
tion to the oo-prolonged PDE system too : £°° C J°°F- Thus, the differential 
forms Q*{£°°) can also be given the structure of a bi-differential bi-graded al- 
gebra (f7*'*(i^),d^,d^), where n*/{F) ^ n*{£°°). The cohomology groups 
77^'*(dh) are called the characteristic cohomology groups of the PDE system £. 
Characteristic cohomology classes can be identified with important geometric 
and algebraic properties of £- For instance, a representative of an element of 
i?^~^'°(dh) is a non-trivial conservation law of £ [71 flOl [23]. 
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3 Covariant presymplectic current 



With notation following Section [2l a local action functional of order fc on _F is 
a function S[(f>] of sections £ r(i^), 

= / (/■<^)*/:, (6) 

where the Lagrangian density C € n"''^(F) is a purely horizontal n-form on 
j'^F. It is local because, given a section </> and local coordinates {x'^,Uj) on 
j'^F, the pullback at x G M can be written as 

{{fcj,rc){x)^c{x\dir{^)), (7) 

which depends only on x and on the derivatives of at x up to order k. The 
integral over M will be considered formal, since all the necessary properties will 
be derived from C. Incidentally the usual variational derivative of variational 
calculus can be put into direct correspondence with the vertical differential dv 
on this complex, which is how the name variational bicomplex was established. 

In Section [5] we introduced the variational bicomplex {Q'^''"{F),dh,dv) of 
vertically and horizontally graded differential forms on J°°F. A Lagrangian 
density, being of top horizontal degree, is then a closed element of ri"'"(F). 
Usually, Lagrangian densities are considered equivalent if they differ by a hori- 
zontally exact term, for example C and C + d^B, where d^B is often referred to 
as a boundary term. In other words, we should think of C not just as an element 
of r2"'°(F), but rather a representative of a cohomology class [C] S i/"'°(dh). 

Note that even though C can be thought of as a form on J'^F, below we 
will carry out all calculations on J°°F, with the proviso that all intermedi- 
ate formulas could have been projected onto jet bundles of some finite order, 
bounded throughout the calculation. When necessary, we shall make use of a 
local coordinate system (x*,uf) on J°°F. 

Using integration by parts if necessary, we can always write the first vertical 
variation of the Lagrangian density as 

dv/: = ELaAdvu"-dh6'. (8) 

All terms proportional to dvwj, |/| > 0, have been absorbed into d^O- In the 
course of the performing the integrations by parts, the coefficients ELa can 
acquire dependence on jets up to order 2k. The form 9 is not uniquely specified, 
as one can freely substitute 9^9 + dhU, where a is any forrrt^ in r2"~^'^(F). 
Thus, the jet order of 9 is not bounded from above (due to the arbitrariness in 
a). However, as can be seen from integration by parts in local coordinates, 9 
can always be chosen to depend on jets of order no higher than 2k — 1. The 
forms ELa define a bundle morphism 

EL : J^'^F ^ F* = A"M ®m F* , (9) 

^Actually, d(,(T could also be replaced by any merely closed form. However, the cohomology 
groups H*'^(F) for v > are always trivial 3, Ch.5], so there is no loss in generality. 
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where F* ^> M is the dual vector bundle to F ^ M . The equation form 
(EL,F*) determines the so-called Euler-Lagrange (EL) PDE system l: ^el C 
J^*'F associated with the Lagrangian density C or equivalently the local action 
functional S'[0]. 

A PDE system with an equation form given by Euler-Lagrange equations 
of a Lagrangian density is said to be variational. The form 9 G n'^~'^'^{F) is 
referred to as the presymplectic potential current density. Applying the vertical 
differential to 9 we obtain the presymplectic current density (or the presymplectic 
current density defined by C if the extra precision is necessary) . 

w ^ dv6l, (10) 

with w e n"^^''^{F). Since the vertical differential does not increase the jet 
order, the jet order of w is bounded by that of 9. Note that changing £ by a 
boundary term d^B changes 9 by the vertically exact term — dyB. Hence, w is not 
altered by this change. If one changes 9 by the addition of a horizontally exact 
term dh cr, then the presymplectic current density also changes by a horizontally 
exact term, w ^ w — dhdvcr. 

The terminology for uj comes from classical field theory in physics. Restricted 
to solution of £el and integrated over a closed codim-1 surface in M (most often 
a Cauchy surface) uj defines a presymplectic form on the space of solutions of 
the EL equations. This (possibly infinite dimensional) space, known as the 
phase space, is then a presymplectic manifold. If the initial value problem on 
this surface of integration is well posed, it is even a symplectic manifold. This 
method of construction the symplectic form on the phase space of a classical 
field theory is known as the covariant phase space method |25[ IHl 1111 119[ |5] . 

Since w is, like A, to be integrated over a boundaryless submanifold, we 
only care about its equivalence class [w] G fl"~^'^ (F) / dhfl"'^^'^ (F) . Thus, the 
discussion above has shown the following 

Proposition 1. Given a Lagrangian density C G il"''^{F) of order k, there exist 
forms 9 £ 51"^^'^(i^) and uj G il"~^''^(_F), of jet order at most 2k — 1, as well as 
a section EL : J^*^ — >■ F* that satisfy the equations 

d^C = ELa A dvu? - dh9 (11) 
d^9 = u. (12) 

Moreover, the equivalence class [C] G -//"'"(dh) determines the section EL and 
the equivalence class [lo] G f2"^^'^(F)/dhf2"^^'^(F) uniquely. 

The following lemma is an easy consequence of the definition of lo. 

Lemma 2. When pulled back along the inclusion l: ^el ^ J^^F, the image 
L*oj G r2^~J'^(i^) of the form LJ Equation ()10|) is both horizontally and vertically 
closed: 

d^EL^*^^0, (13) 
d^'=n*cj = 0. (14) 
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Proof. The horizontal and vertical differentials on iPel are defined by puUback 
along t, that is, df^^^c* = i*dh and d^^^^i* = i*dv. Since w = d^9 is already 
vertically closed in ft'^~^''^{F), it is a fortiori vertically closed in rig~^^'^(F). 
The rest is a consequence of the nilpotence and anti-commutativity of d^ and 
d,: 

= dlC = dvELa A dvu° - dvdh6', (15) 
dhcj = dhdv6l = -d^du9 = -d^ELa A d^u", (16) 
df'=^i*a; = i*dhw = -t*dvELa A dvit° = 0, (17) 

where the last equality holds because ELa and dyELa generate the ideal in 
annihilated by the pullback l*. □ 

In fact, we will promote the name presymplectic current density to any form 
satisfying these properties. 

Definition 1. Given a PDE system l: £ C J^F we call a form u G VL'^^^''^ [F) 
a presymplectic current density compatible with £ if it is both horizontally and 
vertically closed, 

diuj = 0, (18) 
d^Lu = 0. (19) 

In other words, u represents a cocycle (that is, a vertically closed element) in 
the cohomology complex (i/g~^'*(dh), dv). (Note that we are not bounding the 
jet order of tj.) 

The particular form defined by Eq. (jlOp will be referred to as the presym- 
plectic current density associated to or obtained from the Lagrangian density 
C, if there is any potential confusion. 



4 Inverse problem 

In the preceding section we have defined variational PDE systems showed that 
each one is endowed with a geometric structure (the presymplectic current). 
The inverse problem of the calculus of variations (or the inverse problem for 
short) is, given a bundle F — > M and a PDE system l: £ C J^F of order k, to 
decide when it is variational. 

The simpler non- multiplier version of the inverse problem presupposes that 
we are given an equation form {f,F*) for £. It consists of deciding whether 
there exists a Lagrangian density whose Euler-Lagrange equations (EL, F*) are 
equal to {f,F*). The necessary and sufficient conditions for the non- multiplier 
inverse problem are known and are called the Helmholtz conditions. They can 
be formulated elegantly as the requirement that the form fa A dvM° G ri"'^(i^) 
be closed in a slightly extended version of the variational bicomplex [3l Ch.5], 
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The harder multiplier inverse problem consists of deciding variationality di- 
rectly from the sub-bundle l: E C J^F itself or, equivalently, any regular equa- 
tion form (/, E) of £. The name comes from the possibility of reducing it to the 
simpler problem by finding the right set of "multipliers" e (which could also be 
differential operators) such that (e o f,F*) satisfies the Helmholtz conditions. 
Unfortunately, the multiplier inverse problem does not yet have a satisfactory 
solution in full generality [U [21] . 

An important contribution to the subject was made in |15| . as discussed in 
the llntroductionl Henneaux showed that, for second order ODE systems that 
can be put into the canonical form 



the existence of a symplectic form Lj(t,q,q), defined on the bundle of initial 
data {q, q) over the time axis, that is conserved by the flow of the vector field 
associated to the ODE system (^0]). 



is equivalent to this ODE system being variational with a unique Lagrangian 
density C (up to addition of boundary terms) whose associated symplectic form 
is equal to w. Since the EL equations of C will in general not be directly in the 
canonical form (I^U)) . this result shows that a conserved symplectic form w is a 
certificate of a positive solution of the multiplier inverse problem. Henneaux's 
proof even provides a procedure to construct £ from uj and the ODE system. The 
multiplier inverse problem, is then reduced to identifying conserved symplectic 
forms, which could be attacked by algebraic means. 

Unfortunately, until rather recently, it has not been clear how to generalize 
Henneaux's reformulation of the multiplier inverse problem to PDEs. Several 
aspects of the discussion in the previous paragraph are specific to ODEs: (a) the 
possibility of a simple canonical form like (I^Hl) . (b) the geometric formulation 
of the ODE as a vector field, (c) a local symplectic form ui, (d) the conser- 
vation condition (|2ip. In this section, we present a partial generalization of 
Henneaux's result to PDE systems. The PDE analogs of the key aspects are 
(b) the geometric formulation in terms of jet bundles (as in Section [2]), (c) the 
local covariant presymplectic current density lu (as in Definition [1]), and (d) the 
closure condition d^cj — 0. Unfortunately, we have not been able to identify 
simple, local analogs of the canonical form ([20]) and nondegeneracy of ui (that 
is, being symplectic rather than just presymplectic). Due to the last caveat, the 
procedure given below does not produce a unique class of equivalent Lagrangian 
densities C associated to a given PDE system £ and a presymplectic current uj. 
On the other hand, each Lagrangian density produced is in a certain sense a 
subsystem of £: any solution of £ also solves the corresponding EL equations. 
Section [5] is an attempt to characterize the class of Lagrangian densities that 
can be so produced. 

The key observation that connects a local presymplectic current density 
with a variational formulation was made in fO', Sec. 4], which is a more geometric 



Q - f{t,q,q) = 0, 



(20) 



dtUJ — J>ffU] = 0, 



(21) 
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formulation of an earlier observation made in [TB]. However, these authors did 
not attempt to place this result in the context of other work on the inverse 
problem of the calculus of variations and did not remark the similarity with 
the previous work of Henneaux. Moreover, their calculations remained "on- 
shell", which avoided lifting the Lagrangian density "off-shell" (see the proof 
below), which is really necessary for a solution of the inverse problem. Below, 
we clarify this observation and show in detail how a Lagrangian density can 
be constructed using a method related to cohomological descent (see also 
Refs. [88,89] and [191] therein). 

As before, consider a vector bundle F ^ M over an n-dimensional manifold 
M and a regular PDE system t: £ C J'^F oi order k. Recall also that F* = 
A^M (g)M F* is the densitized dual vector bundle of F. Finally, an important 
hypothesis currently assumed is that the de Rham cohomology groups of £°° 
and J°°F all vanish. All other relevant notions and notation are defined in 
Section m 

Theorem 3. // there exists a presymplectic current density oj G fl^ ^'^{F) 
compatible with £, then there exists a local Lagrangian density C € ri"'*'(_F) 
such that the associated presymplectic current density u) coincides with u) on 
solutions (i*^ui — p*ao'^) o.i^d, the Euler- Lagrange PDE system £'el, with equation 
form (EL, _F*), is compatible with £ (£°° C £^^ or all solutions of £ also solve 
£el)- 

Proof. The fact that uj, from the definition of a presymplectic current density, is 
both horizontally and vertically closed as an element of rj^"^'^(i^) can be easily 
seen to be equivalent to lj being de Rham closed as an element of f2g^^(F). 
Equivalently, if w is of homogeneous degrees (n — 1, 2) and de Rham closed, it 
follows that it is both horizontally and vertically closed. 

By assumption, the de Rham cohomology group H^^^ (F) is trivial. If it were 
not, there could be global topological obstructions to this construction, which 
we do not discuss in the current treatment. Then there must exist an element 
p G ^lg{F) such that d^/5 — Gj. If we expand p in components of homogeneous 
vertical and horizontal degrees, we can represent it as p = y^^^^,_^ Ph.v This 
sum is finite, since < /i < n and < The equation 

d^p = d^p + d^/5 = c2; (22) 

then naturally expands into the following system for the homogeneous compo- 
nents 

/5o,„ = 0, (23) 
dv/5i,„-i = -d^/5o,„, (24) 

: (25) 

dv /0«-i,i = -d^/5„_2,2 + (26) 
d^/5„,o = -d^p„_i,i. (27) 
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Note that p could be constructed by solving the above equations term by term. 
This method is a special case of the cohomological descent method [7]. 

As written, these equations hold "on-shell," that is, on the oo-prolonged 
PDE manifold £°°. But the descent equations lift "off-shell," to the total space 
of jet bundle J°°F containing £°° . Denote the lifted forms by removing hats, 
Ph,v = t'laPh,v, with the exception & — t^w'. These lifts are not unique, as we 
could always change them by adding terms that are annihilated by the puUback 
to £°°. RecaU that £°° is defined by the equations p°°f = 0, so given local 
coordinates {x\va) on the equation bundle E, which extend to {x\via) on 
J°°E, the terms annihilated by the pullback to £°° must be proportional to fiA 
or the exterior vertical derivatives dv//A- After the lift, the above equations for 
Ph,v also only hold up to terms proportional to Jja or dv//A, 

dvPO,n = + flA^o^i+l + dvfiA A /Ig^j, 

dvPl,n-l = — dhPo.n + flA^i^ + dvfiA A ^{^^1, 



dvPn-2,2 = -dhPn-3,3 + flA>^i^2,3 + dv//A A pi^2,2^ 

and 

dvPn-ia = — dhP„-2,2 + flA^i^ 1^2 + dy/zA A l-J-i^iA + "^'i 
dvPn.O = ^dhP„-l,l + flA^n^l + ^vflA A fin^fi- 

Note that we have introduced the coefficient forms A^"^ and fi^f^^ to parametrize 
the terms annihilated by the pullback to £°° ; they are of homogeneous horizontal 
and vertical degrees, as indicate by their subscripts. These coefficient forms are 
not simply arbitrary. As shown below, Xj/^i and will contain in them the 
information about the multipliers needed to solve the inverse problem. 

Given local coordinates {x^,uf) on J°°F, the goal now is to construct the 
forms £, 9, uj and EL^ A dvM°, of respective degrees (n, 0), [n — 1, 1), (n — 1, 2) 
and (n, 1), such that the corresponding equations of the covariant phase space 
method hold: 

dv/: = ELa A dvu° - dh6i, (28) 
dv6l = UJ. (29) 

We rewrite the (n, 0)-descent equation as 

dv(Pn,0 - flApij)) = flAiX'rA " dvM^^o) " dhP«-l,l (30) 
= flAe'/ A dvU" - du{pn-lA + flAX'n-l.l). (31) 

where integration by parts was used to construct X'n-i,! and e^^, with the latter 
being (n, 0)-forms. We add dv(//AA^:5i i) to both sides of the {n — 1, l)-descent 
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equation and rewrite it as 



dv(Pn-l,l + flA^n-l,l) = - dhP„-2,2 + //a(A^"1i 2 + '^v^n^l^l) 

+ dJjAA{fii\,+X'^\,). (32) 
It is now clear that these equations take the desired form with the following 



identifications: 

^ = Pn,0 " //AAifiJ), (33) 

O^Pn-lA+flAX'^^l,l, (34) 

UJ = uj' - dhp„-2.2 (35) 

ELaAd.u'' = fiAe'/d.u''. (36) 



The form ELa A dvM° naturally correspond to a bundle morphism EL: J F ^ 
F*, for some finite jet degree /, which defines the Euler-Lagrange PDE system 
Eel C J''F via the equation form (EL, F*). From the last equation, it is obvious 
that the constructed Euler-Lagrange PDE system contains the original one, 
C £^^. □ 

Note that all the above calculations were done on jet bundles of infinite 
order. However, each of the forms introduced at intermediate steps depends 
only on jet coordinates of some finite order. Therefore, the final Lagrangian 
will also depend on jet coordinates up to some finite order, which may be much 
higher than the order of the original PDE system. Note that the degree of the 
Lagrangian could be artificially infiated by the presence of boundary terms like 
dhS, where B could depend on jet coordinates of some high order. 

It is also important to remark that the resulting Euler-Lagrange equations 
may not be equivalent to the full original PDE system, but only to a subsystem 
thereof or a "weaker" system, one whose solution space contains all solutions 
of £, but may be strictly larger. This is unavoidable, since the PDE system 
may consist, for example, of several uncoupled subsystems, each of which may 
have an independent variational formulation. More complicated situations are 
of course possible. In the ODE context, the requirements that uj actually be 
symplectic (rather than just presymplectic) and that the ODE system is in the 
canonical form are sufficient to guarantee that the original system is fully 
variational [T5l. Unfortunately, at least one of these conditions fails already 
when the ODE system is not determined (possibly because of gauge invariance) 
or under the inclusion of algebraic equations that constrain the initial data [q, q). 
In the PDE case, in analogy with the ODE one, there may be a set of conditions 
on the symbol of the PDE system and on uj that guarantees a fully variational 
formulation, which may be further complicated by allowing equations with gauge 
symmetries or constraints. However, such analogous conditions remain to be 
investigated in detail. 
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5 Arbitrariness in Lagrangian density construc- 
tion 

There were a number of choices involved in the construction of the (off-shell) 
Lagrangian density H from the (on-shcU) prcsymplcctic current density a). In 
this section, we investigate how the resulting £ depends on these choices. 
The choices are exhausted by the following substitutions: 

(i) Cj^Cj^ d^7r„_2,2, with d57r„_2,2 = 0; (37) 

(ii) p^.p + d^CT; (38) 

(iii) p^p + Z/A^^+dv/zAA/x^^, (39) 
a;'^a;' + /jAA'^^+d,/7AA/x'^^. (40) 

We deal with each kind of substitution one by one. 

(i) It is easy to see the following identity: 

d^(/3 + ^„-2,2) = w + d^7r„_2,2. (41) 

Since the change p — p-|-7r„_2,2 affects neither of the p„,.o or Pn-i.i components, 
the off-shell lift of the (n, 0) descent equation is unmodified. Therefore, the 
Lagrangian density L does not change. 

(ii) Lifting u off-shell to u and using the decomposition into homogeneous 
components, u = "Yiih+v^n '^h,v, we see the change 

Ph,v — >■ Ph,v + dhC/i-l.u + dvCT^t,„_i. (42) 

However, since d(p + dcr) = dp, the descent eqiiations arc; unmodified. In the 
end, the Lagrangian changes only as £ ^ £ + dhCT„_i,o- Since the change is by 
a horizontally exact term, the EL equations remain the same. 

(iii) This kind of substitution does in general change the equivalence class 
of the Lagrangian density, that is, the EL equations of the modified Lagrangian 
density may be different. The forms parametrizing the failure of the off-shell 
lift of the descent equations undergo the change (except for the Xn-1,2 and 
Hn-1,1 coefficients, which undergo an obvious additional changes compensating 
the change in ui', which ultimately does not affect £) 

^h.^ ^ + (dh Aj^fi,. + dx^^ A X'JXv) + dvAl^li , (43) 
<^ ^ <^ - (dhpi/-^i.„ + dx(^ A f^i^Xv) + C^h,v - dvpj,^ti). (44) 
Therefore, the Lagrangian density undergoes the change 

£^£ + £ = £ + fuA~Xi!,^ - fiiA [At^^ - (dhpJ.^A.o + dx^' A fi'^Xo)] (45) 

= £ + fiiA{d^,llt!A,o + dx^' A (46) 

In general, the multipliers e^^ will change as well, say e^'^ — ?■ e^'^ -|- e^'^. Thus, 
the EL equations of £ and £ + £ may not be equivalent. In other words, the 
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equivalence classes [C] and [C + C] will differ. Though both sets of EL equations 
will be consequences of £. 

At this point, it is worth reflecting on when two Lagrangians C and C should 
be considered equivalent. The standard answer is iff they differ by a boundary 
term, C — C = d^B. However, consider the simple 1-dimensional Lagrangians 



where a is a constant (though in principle it could be a more complicated 
function of A and qi). Note that the difference between the two Lagrangians, 
£' — £- = a{qi — 92), is not a boundary term, though it is proportional to the 
constraint equation (72 ~ 91 = obtained by varying A. It is easy to check that 
the resulting Euler-Lagrange equations for either Lagrangian are equivalent to 
the set qi = 0, q2 = qi, X ~ 0. Moreover, their symplectic currents agree as 
well, uj — dv(?i A dvQi = uj' . Furthermore, since we never limited the order 
of the Lagrangian density, it could be of order zero (a ordinary, non-diffcrcntial 
variational problem). Then both w = a;' = and the critical points of C and C 
coincide as long as dv(-C' — C) vanishes at all critical points. Therefore, in the 
presence of constraints on the initial data, the classes of equivalent Lagrangian 
densities are larger than just those that differ by boundary terms, if equivalence 
is evaluated in terms of the Euler-Lagrange equations and the presymplectic 
current. 

Consider the following preordeiH on Lagrangian densities. We say that C 
contains £, £ ^ if the 00-jet bundle submanifolds defined by the prolonga- 
tions of the corresponding Euler-Lagrange equations obey £^ C and their 
presymplectic current densities agree up to horizontally exact term when pulled 
back to the more restrictive equation, lI^lj = l^^uj' , where l: £el C J°°F. The 
relation -< is clearly transitive and reflexive and hence a preorder. 

The containment relation is however not a partial order, because it is not 
antisymmetric. As already discussed earlier, ii C — C — d^^B, then L < C 
and L' -< C, while C ^ C . On the other hand, this preorder gives rise to the 
equivalence relation C, ^ CJ (iff £ ^ C and CJ < C) and projects to an actual 
partial order [C] -< [C] on the equivalence classes with respect to ■~. 

Given a PDE system t: £ C J''F with a compatible presymplectic current 
density uj, define the set Lg^ci to consist of equivalence classes of Lagrangian 
densities [C] such that £°° C £^ and l'^uj — p*^oj. Then Ls^q is an upper set 
with respect to the -< partial order (which means that [C] £ Lg.cD and [C] -< [C] 
implies [C] G Ls^cj). After all, as should be clear from the definition of the -< 

preorder P on a set X is a relation such that is reflexive (xPx) and transitive (xPy 
and yPz implies xPz), but in general neither symmetric nor antisymmetric. All partial orders 
and equivalence relations are preorders. The maximum symmetric subrelation E {xEy iff xPy 
and yPx) is an equivalence relation and the quotient X — X/ E projects P to a partial order 
P/E on X/E (which is necessarily antisymmetric) I22| . 



'^^['7i,'Z2, A] = T^ql + ^ql + \{qi - qi) dt, 



(47) 



C'[quq2,M^ 2*^1 + 2'^^ + (^ + ") (91 -92) dt, 



(48) 
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relation, if p = C + 9 could arise in the construction in the preceding section 
(that is £ G L£_c), then so could p = C + 9' (that is C £ Ls^q) for any £' that 
contains £{£-<£'). 

It is clear that the construction of the preceding section can yield any ele- 
ment of as a partial solution of the inverse problem. It is possible that a 
consideration of the structure of the partial order -< as well as the decompos- 
ability and minimal elements of Lg^^ in relation with the symbol of the PDE 
system £°° can yield a definite solution of the inverse problem. These questions 
are yet to be investigated in detail. 

6 Conclusion 

We have shown, following a recent observation in [SI [TB] and a strong analogy 
with previous work in [1^1 [HI [T] , that a horizontally conserved presymplectic 
current density is a certificate that a subsystem of a given PDE system or a 
comparatively weaker system is variational (any solution of the original PDE 
system also solves the obtained variational system) . If this subsystem is actually 
the full system, then the multiplier inverse problem of the calculus of variations 
has a positive solution. 

Restricting to the context of second order ODEs that can be put in canon- 
ical form ([^O)) . it is known that any two Lagrangians whose Euler-Lagrange 
equations are equivalent to the given ODE system £ and that have equivalent 
symplectic current densities lu must differ by a boundary term[15]. In that case, 
the equivalence class of Lagrangians solving the full inverse problem of the cal- 
culus of variations for £ depends only on the characteristic cohomology class 
in [lu] e i7^'~^'^(dh). However, in more general contexts (for general PDEs, 
or even ODEs with constraints) a characterization of equivalent Lagrangians 
(those sharing the same set of solutions and equivalent presymplectic current 
densities) is still missing. 

On the other hand, in Section [Sj we defined a preorder relation -< on La- 
grangian densities C G f2"'°(F) that relates both this equivalence problem for 
Lagrangian densities and the ambiguity in the Lagrangians produced by the 
construction of Section |4| The precise structure of the preorder ~< and the con- 
ditions on a PDE system and compatible presymplectic current that guarantee 
that it is fully variational (rather than just a subsystem thereof) remain to be 
investigated. It is possible that generic methods for determining the characteris- 
tic cohomology groups of a PDE [TUl [23] would be helpful in classifying possible 
presymplectic current densities and hence variational formulations. 
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